Vertical loop nodes in iron-based superconductors 
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We consider Fe-based superconductors with s H gap with accidental nodes on electron pockets. 
We analyze how the gap structure changes if we include into the consideration the hybridization 
between the two electron pockets (the inter-pocket hopping term with momentum (7r,7r,7r). We 
derive the hybridization term and relate it to the absence of inversion symmetry in the Fe-plane 
because of two non-equivalent locations of pnictogen (chalcogen) above and below the plane. We 
find that the hybridization tends to eliminate the nodes - as it increases, the pairs of neighboring 
nodes approach each other, merge and disappear once the hybridization exceeds a certain threshold. 
The nodes disappear first around k z = n/2, and vertical line nodes split into two vertical loops 
centered at k z = and k z — n. We also show that the hybridization moves the nodes along 
the loops away from the normal state Fermi surfaces. This creates a subset of k— points at which 
the peak in the spectral function does not shift as the system enters into a superconducting state 
("no-shift" lines). These "no-shift" lines evolve with increasing hybridization in highly non-trivial 
manner and eventually form horizontal loops in (k x , k y ) plane, surrounding the nodes. Both vertical 
line nodes and horizontal "no-shift" loops surrounding them should be detectable in photoemission 
experiments. 



I. INTRODUCTION 

Understanding high-temperature superconductivity in 
doped Fe-pnictides and Fe-chalcogenides remains the top 
priority for the condensed matter community pJ-BJ. Su- 
perconductivity in weakly /moderately doped systems is 
generally believed to be the consequence of the complex 
geometry of the Fermi surface (FS), which consists of hole 
and electron pockets located in different regions of the 
Brillouin zone. The prevailing scenario is that the super- 
conducting gap has an s— wave symmetry, but changes 
sign between hole and electron pockets [5, 6], and may 
even have accidental nodes [E0]- 

Previous studies of the pairing mechanism and the gap 
structure in Fe-based superconductors mostly focused on 
a 2D model of a single Fe plane with adjacent pnicto- 
gen/chalcogen atoms located immediately above and be- 
low it [ll.l61.l8l4l3j. Electrons from 3d-orbitals of Fe hop 
mostly indirectly, via 4p pnictogen/chalcogen sites. The 
low-energy structure around Ep obtained from a 3d-4p 
hopping has been fitted fll. Ifl IR [TlTjl3| to an effective 
"Fe only" 2D 5-orbital tight-binding Hamiltonian with 
inter-orbital and intra-orbital hopping terms. The band 
structure and the location of the electron and hole FSs 
have been deduced by evaluating the eigenvalues of the 5- 
orbital tight-binding Hamiltonian and analyzing the en- 
ergy profile. The pairing problem can be most straight- 
forwardly analyzed in the band basis, by associating the 
operators corresponding to eigenfunctions with band op- 
erators and re-writing the interaction in the band ba- 
sis The band operators are linear combinations of 
orbital operators, and the interactions in the band basis 
are the ones in the orbital basis, dressed by "coherence 
factors" associated with the fact that each eigenfunction 
is a linear combinations of Fe-orbitals. The dressed inter- 
actions inherit angular dependencies from the coherence 
factors. Solving for gaps on different FSs, one then gener- 



ally obtains angular-dependent gaps, even in the s-wave 
case. By generic reasoning, the s— wave gaps on hole 
pockets contain cos4rt</> harmonics (n = 0,1,2..) where 
4> is the angle along the hole pocket, while the gaps on 
electron pockets contain both cos An8 and cos26*(2n + 1) 
harmonics, with the angle 8 counted relative to, e. g., x 
axis for one pocket and y axis for the other 
Numerical analysis shows [JUII that cos 29 harmonic of 
the gap is the strongest one, and in some materials it is 
large enough to induce accidental nodes on electron FSs. 

The actual BZ, however, contains two Fe atoms be- 
cause pnictogen/chalcogen is located below and above 
Fe plane in checkerboard order. The positions of pnicto- 
gen/chalcogen atoms are shifted by half of lattice spac- 
ing in both x and y directions relative to the positions 
of Fe. A half of pnictogen/chalcogen atoms are located 
above Fe plane and half are located below Fe plane. As 
a consequence, for half of Fe atoms hoping to pnicto- 
gen/chalcogen means hopping down along z and for the 
half it means hopping up along z. Therefore, the sym- 
metry of the original lattice is lower than the symmetry 
of the Fe lattice, and the effective tight-binding Hamil- 
tonian should generally contain two type of terms - the 
ones with zero momentum transfer, and the ones with 
momentum transfer Q — (ir, it) [l6l . [l7| (here and below 
we set interatomic spacing to one). 

The tight-binding Hamiltonian considered in earlier 
studies [l| includes only the terms with zero momen- 
tum transfer. For such a Hamiltonian, the transforma- 
tion from 1-Fe BZ to 2-Fe BZ is just a rotation in a mo- 
mentum space: momenta k x and k y from 1-Fe zone are 
transformed into k a — k x + k y and k^ — k x — k y in the 
2-Fe zone. The gap structure is obviously not affected 
by this transformation, up to a change of variables. In 
particular, if the gap has accidental nodes in 1-Fe zone, 
it will have accidental nodes in 2-Fe zone as well. 

In this paper we demonstrate that the effective tight- 
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FIG. 1: (a) Schematic representation of a single layer of pnic- 
tide/chalcogenide. The solid square designates a Fe (solid cir- 
cles) only lattice unit cell. The As/Se atoms form two inequiv- 
alent subsets. The atoms denoted by a cross and empty cir- 
cles are shifted off the Fe plane in opposite directions, (b) The 
hopping amplitudes to a neighboring pnictogen/chalcogen site 
are generally different for pnictigen/chalcogen above and be- 
low Fe plane, i.e., in general A\ 7^ A2. 



binding Hamiltonian also contains hopping terms with 
momentum transfer Q, and that these terms affect the 
physics in a quantitative way, and, in particular, sub- 
stantially modify the gap structure. The (tt, tt) term 
vanishes if one approximates a 2D 5-orbital model by 
just d xz and d yz orbitals, as has been done in some ear- 
lier studies dUEII but is present once one adds a d xy 
orbital (the (n, 7r) term describes the hopping from d xy 
orbital to d xz or d yz orbital). Because the low-energy 
hole and electron states in the full 2D tight-binding 
Hamiltonian have contributions from d XZl d yz , and d xy 
orbitals [ljjg, HffldEGl, the (tt.tt) hopping term sur- 
vives the transformation to band formalism and becomes 
iptjpk+Q) where ip are band operators. As the two elec- 
tron FSs are separated exactly by Q, such term gives rise 
to a hybridization between the two electron pockets. The 
new band operators for low-energy fermions then become 
linear combination of the original "0fc , and the pairing in- 
teraction, re-expressed in terms of the new band opera- 
tors, acquires additional coherence factors, which modify 
its angular dependence. The modification of the angular 
dependence of the interaction in turn modifies the angu- 
lar dependence of the gaps on electron FSs. This physics 
must be present in 10-orbital studies d, Qii [H| , which 
treat orbital excitations with momentum k and fc + Q as 
separate degrees of freedom, but we are not aware of the 
attempt to fit the results of 10-orbital numerical studies 
by an effective tight-binding Hamiltonian with an addi- 
tional (7r, 7r) term. 

In the previous work [2pj| we demonstrated that hy- 
bridization plays a crucial role in systems with only elec- 
tron pockets, like KFe2Se2. Several groups analyzed 
the pairing problem in KFe2Se2 in the absence of hy- 
bridization. If the interaction between the two elec- 
tron pockets is repulsive, the pairing is possible if inter- 
pocket interaction exceeds intra-pocket repulsion, but the 
gaps on the two electron pockets should have opposite 
signs 0, HH, H3 ■ Such a state has a d— wave symme- 



try because the gap changes sign under rotation by 90°. 
The hybridization mixes the two original electron pockets 
into two new electron pockets, each contains states from 
both original pockets. The two new pockets split upon 
hybridization, and the solution of the pairing problem 
for strong enough hybridization shows that the system 
prefers to form an s— wave gap with opposite sign on the 
two hybridized FSs |17ll20j . In terms of original fermions 
the condensate wave function for such an s-wave state 
is made out of fermions from different electron FSs, i.e., 
a pair has a momentum Q. In terms of new, hybridized 
fermions, such a state is a conventional one, with zero to- 
tal momentum of a pair. The microscopic analysis of the 
pairing in the presence of hybridization shows [2(| that, 
with increasing hybridization, the system undergoes a 
two-step transition from a d— wave state to a d±is state 
and then to an s" 1 state. 

In this paper we analyze the role of hybridization in 
systems in which both hole and electron pockets are 
present and the gap has s H form. We take the gap 
structure from the 2D tight-binding Hamiltonian without 
Q terms as an input and analyze how it changes once we 
add additional hopping with momentum transfer Q. At 
first glance, the hybridization between the two electron 
pockets should not lead to qualitative changes in the s" 1 
gap because the gaps order parameter on these pockets 
on average have the same sign. We show, however, that 
the actual situation is more tricky and hybridization does 
give rise to qualitative changes in the gap structure in 
systems like BaFe^Asi-^P^ (Ref. @), LaOFeP (Ref. 
[23 |) and LiFeP (Ref. (25|), in which s H gap has acci- 
dental nodes, which most likely reside on electron pock- 
ets. By symmetry, there are four nodes on each electron 
pocket (8 nodes in total) . We show that the hybridization 
brings the neighboring nodes close to each other. When 
the hybridization reaches a certain threshold, the nodes 
coalesce and disappear (see Figs. I2l4p . 

The importance of hybridization between electron 
pockets was recognized in Ref. The authors of that 
work conjectured that it may affect the location of the 
gap nodes. Our explicit calculations do show that the 
nodal structure is indeed affected by hybridization. 

We follow the evolution of the nodes before the hy- 
bridization reaches a threshold. We show that at any 
non-zero hybridization the nodes do not reside on the 
normal state FSs but are shifted into the fc— region be- 
tween the hybridized FSs, where the ARPES intensity 
above T c is peaked at a finite negative frequency. As 
a result, the peak in the ARPES spectra at the nodal 
point doesn't stay intact, as in a "conventional" super- 
conductor with gap nodes, but shifts towards a smaller 
frequency below T c . We further show that each nodal 
point is surrounded by "no-shift" lines (NSL), at which 
quasiparticle energies in the normal and superconduct- 
ing states are equal, i.e., the peak in the ARPES spec- 
trum does not shift when the system becomes a super- 
conductor (Figs. Without hybridization, these NSL 
are radial beams transverse to the FS, in the direction 



3 



along which the gap vanishes. Upon hybridization, NSL 
rapidly evolve, rotate by 90°, and transform into loops 
directed along the FSs. Inside the loop, ARPES peak 
moves to a lower negative frequency below T c , outside 
the loop it moves to a higher negative frequency. We pro- 
pose to search for such NSL in ARPES measurements on 
BaFe2(Asi_ x P x )2 and other s H superconductors with 
gap nodes, like LaOFeP ^ and LiFeP [Hj]. 

We next consider the effect of hybridization in 3D sys- 
tems. For 1111 systems with a simple tetragonal lattice, 
the hybridization term has 3D momentum (n, tt, 0), 3D 
effects are secondary, and the consideration is qualita- 
tively the same as in 2D systems. In 122 systems with 
body-centered tetragonal lattice, the situation is more in- 
teresting. We show that in these systems the hybridiza- 
tion between electron pockets is three-dimensional, with 
hybridization vector (7r, 7r,7r), even if we neglect the k z 
dependence of the gap along the third direction within 
the lFe tight-binding Hamiltonian (26[. This is due to 
the fact that in 122 materials pnictogen/chalcogen atoms 
alternate in checkerboard order between neighboring Fe 
planes, i.e., if for a given Fe plane a pnictogen/chalcogen 
atom at a given (x, y) is located above Fe plane, then 
in the neighboring Fe plane a pnictogen/chalcogen atom 
at the same (x,y) is located below Fe plane [271 ]. We 
show that, for such a structure, the hybridization term 
has 3D momentum Q = (7r,7r, tt), and is the largest at 
k z = tt/2 (it actually does not vanish, except for k z — 
and k z = tt, even if we consider only d xz and d yz Fe 
orbitals). We show that, at the same time, the critical 
value of the hybridization, at which the nodes disappear, 
is the smallest at k z = tt/2 and is much larger for k z = 
and k z — tt. The combination of the facts that at tt/2 the 
value of the hybridization is the highest and the thresh- 
old value is the lowest implies that there must be a wide 
range of hybridizations when the nodes are eliminated 
near tt/2 but are still present near k z = and k z = tt. 
In this range vertical line nodes close up into vertical 
loop nodes centered at k z = and k z — tt (see Fig. [5]). 
The gap structure with such loop nodes has been pro- 
posed phenomenologically [23j as the best candidate to 
fit the thermal conductivity data in BaFe2(Asi_ a: P :r )2, 
for which penetration depth, thermal conductivity, spe- 
cific heat and NMR data all show [28[ that the gap must 
have nodes, and at least some ARPES results [29| indi- 
cate that the nodes must be on the electron FSs. Our 
result provides a microscopic explanation of vertical loop 
nodes. 

The paper is organized as follows. In the next section 
we present the model of a s" 1 superconductor with nodes 
on electron pockets and introduce the hybridization am- 
plitude. In Sec. |m] we discuss how the hybridization 
affects the nodes of the superconducting gap in systems 
with a simple tetragonal lattice structure. We consider 
first a toy model with circular electron pockets and then 
consider a more realistic model of elliptical pockets. We 
show that in both cases the neighboring nodes move to- 
wards each other as hybridization amplitude increases. 



In Sec. IIVI we extend the analysis to a 3D body-centered 
tetragonal lattice and show that hybridization splits ver- 
tical line nodes into vertical loop nodes centered at k z — 
and k z — tt. In Sec[V]we discuss the experimental sit- 
uation and comment on recent ARPES studies of the 
superconducting gap structure in BaFe2(Asi_ a ;P x )2. We 
present our conclusions in Sec. I VII We discuss technical 
issues in detail in three Appendices. In Appendix |A1 we 
discuss in detail the microscopic mechanism of hybridiza- 
tion in 2D and 3D lattices and obtain explicit expressions 
for the hybridization amplitude for a simple tetragonal 
lattice structure and for a body-centered tetragonal lat- 
tice structure. In Appendix IBl we discuss the evolution of 
NSL from radial beams in the direction transverse to the 
FS along which the gap vanishes to closed loops directed 
along the FSs. In Appendix [C] we discuss subtleties in 
extracting the positions of gap nodes from ARPES data. 



II. THE MODEL 

We consider an superconductor with two electron 
pockets at (0, tt) and (tt, 0) in the lFeBZ at any given 
k z , and the appropriate number (2 or 3) of hole pockets 
centered at (0,0). The actual pocket structure in Fc- 
based superconductors is somewhat more involved, i.e., 
low-energy states near (0, tt) and (tt, 0) likely contain not 
only electron pockets centered at these points, but also 
hole barrels, centered somewhat away from (0, 7r) and 
(tt,0) (Ref.[30]). These hole barrels will not play a role 
in our consideration and we neglect them. In the absence 
of hybridization, the quadratic part of the Hamiltonian 
for electron pockets is 

H 2 = -I A- ' , A . ' : /. + 6,fc+QV4,fc+Q^2,fc+Q , (1) 
k 

where i/>i2 describe two electron bands, one with low- 
energy excitations near k\ = (0,7r), and another with 
low-energy excitations near k 2 = (tt, 0) — fei + Q (mod- 
ulo 2tt), and ^2,k+Q are the corresponding electron 
dispersions, and Q = (tt,tt). Wc approximate fcrmion 
excitations near the pockets by £i k — vp(9)(\k — k\ \ — 
k F (0)), 6,fe+Q = vf(6 + Tr/2)(\k '- k 2 \ - k F (8 + tt/2)), 
where 9 is the angle along each of the FSs counted from 
the a;— axis [H-[3l.[20T]. By virtue of tetragonal symmetry, 
v F (8) = v F (l + acos29) and k F (9) = k F (l + bcos29). 
The parameter b accounts for the eccentricity (elliptic- 
ity) of the FSs. For 1111 systems, parameters a and b 
are essentially independent on k z , and one can reduce 
the analysis to 2D model. For 122 systems a and b do 
depend on k z and change sign and magnitude between 
k z = and k z = tt. 

We assume, following earlier works [H-flQ, that the 
dominant pairing interaction is between electrons and 
holes, and that its dependence on k z is weak and can 
be neglected. We follow Refs. @,[HL[l5[ and approximate 
electron- hole interaction by a constant term and by cos 29 
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term, which changes sign between the two electron pock- 
ets. Within this approximation, the gaps on hole FSs are 
angle- independent, while the gaps on the two electron 
FSs are A(l ± a cos 26^.). The corresponding effective 
BCS Hamiltonian is 

H BCS [(1 - acos20 k )4>i,k4>i,-k + h.c] 

k 

+ [(1 + a cos 29 k )^ 2 ,k+Q^2,-k-Q + h.c] . (2) 

We consider the case a > 1 , when the gaps have acciden- 
tal nodes in the absence of hybridization between the two 
electron pockets. The hybridization term in the lFeBZ 
has the form 



H Q =y~] Mk)ipl k ip 2 ,k+Q + h.c 



(3) 



We present the derivation of Eq. ([3]) in Appendix [A] for 
both 2D and 3D systems. We show there that the prefac- 
tor A(fc) is generally a complex number, which depends 
on both k z and k±. When only the d xz and d yz or- 
bitals of Fe are considered, the hybridization vanishes 
completely in 2D case and along particular directions in 
3D case. The fc-dependence becomes less strong once one 
adds into consideration d xy orbital along with spin-orbit 
interaction. In the latter case, |A(fc)| does not vanish 
for any k, although it is smaller at k z = and along 
k x = ±k y measured from the center of one of electron 
pockets. The non-singular fc-dependence of A(fc) is not 
essential to our analysis, and to simplify the presentation 
we approximate A(fc) by a constant A for the rest of the 
paper. 

Our goal is to analyze what happens with the nodes, 
and, more generally, with the fermion dispersion, when 
we solve for the pairing in the presence of the hybridiza- 
tion term. 



III. THE EFFECT OF HYBRIDIZATION ON 
THE NODES, 2D CASE 

In this section we consider the case of a simple tetrago- 
nal lattice for which the pairing problem can be analyzed 
within a single 2D cross section (we recall that we treat 
interactions as independent on k z ). In the presence of 
Hq, Eq. ([3])), the quadratic part of the Hamiltonian for 
fermions near the two electron FSs becomes, instead of 

#2=y^ £l,fe/l,fe/l,fe + £,l,k+Qf^ k+Q f2,k+Q 



2,k+Q 



-h.c. 



(4) 



For convenience, we redefine k and count it relative to 
k\ = (0,7r) for /i fermions and relative to k-i = (tt,0) 
for f% fermions, i.e., absorb Q into new k. In these new 



notations, the quadratic Hamiltonian becomes 

j=l,2 k k 



(5) 



To analyze the pairing, we now have to introduce new 
fermions which diagonalize the quadratic form in Eq. (|5|) , 
re-express the pairing interaction in terms of these new 
fermions, and solve for the gaps on electron FSs and 
quasi-particle dispersion. We start with the case of cir- 
cular electron pockets, and then extend the analysis to 
elliptical pockets. 



A. Circular electron pockets 

For circular pockets, = £ 2 ,fc = £,k — VF{k — kp) 
is independent on the angle along the FS. The quadratic 
part of the Hamiltonian is diagonalized by the transfor- 
mation to new operators a and b via /1.2 = (a T b)/^/2 
and becomes 



a\cik 



E 

k 



Zkblbk 



(6) 



with £ fc = ± A. For angle-independent A, the new FSs 
remain concentric circles with different radii. 
The BCS Hamiltonian takes the form 

£/bcs =A ^2 { ( a ka-k + b k b- k ) 
k 

+ {a-kb-k + akb-k)a cos 26 k ) + h.c. (7) 

We see that, in terms of a and b operators, the cos 29 
term measures the strength of inter-band pairing. As A 
increases and the two hybridized FSs become more sep- 
arated, inter-pocket pairing becomes less important, and 
one expects that the angular dependence of the pairing 
interaction will play a lesser role. This is what calcula- 
tions show, as we demonstrate below. 

The quasi-particle dispersion is obtained from zeros of 
the inverse propagator GJ 1 ^ = to I — M with matrix M 
and a unit matrix / operating in Nambu space ipk = 
[ofe, alfe, ^ fc ] tr - We have 



M 



A Ay' 

A -£+ Ay 

Ay £_ A 

Ay A 



(8) 



where y = acos26>fc. The nodal points are located at k r , 
for which det M(k) = 0. Equation ([5]) gives 



det M = 4A 2 e k + $ - A 2 + A 2 (y 2 - 1) 



(9) 



We see that the nodes reside on the "bare" FS, = 0, 
and their angular position 9 k is set by 



cos 9 k = ± 



VP 



oA 



(10) 
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FIG. 2: Color online. The evolution of the FS and gap nodes 
shown by (red) dots for circular pockets and constant (angle 
independent) hybridization amplitude. Without hybridiza- 
tion, the two electron FSs are identical, and each has four 
nodes (panel a). Once the hybridization parameter A becomes 
non-zero, the two electron FSs split, but the nodes remain on 
the original (non-hybridized) FS (panels b and c). As A in- 
creases, the pairs of neighboring nodes come closer to each 
other and eventually merge at A c (panel d) and disappear at 
larger A. We set a — 1.3. 



Equation (| 10[) has eight solutions along a circle (4 pairs 
of nodes near 9 = and other symmetry-related points 
9 = tt/2,9 = 7r, and 9 = 3n/2). As A increases, the two 
nodes located, e.g., near 9 = move toward each other, 
and at a critical A c given by 

A c = AVa 2 - 1 . (11) 

the neighboring nodes merge along symmetry lines. At 
larger A the nodes disappear. We show this behavior in 
Fig. [2j Observe that the critical A c is of order A, i.e., is 
rather small. We will see below that A c is much larger 
when electron FSs are not circular. 

In the presence of impurities, the nodes disappear al- 
ready at A smaller than A c , once the distance between the 
nodes becomes smaller than some minimum value set by 
impurity scattering |3lLl32T ]. 

Consider the range A < A c in more detail. That the 
nodes are not located on the actual (hybridized) FSs has 
a profound effect on the ARPES spectrum. In the nor- 
mal state, the quasiparticle energies are ± A, i.e., one 
energy is positive along = and the other is negative. 
Because ARPES intensity is proportional to the Fermi 
function, ARPES will only detect a negative mode at 
uj = —A. That the nodes in the superconducting state 
are located along = then implies that around these 
fc, the position of the maximum in the ARPES spectra 
shifts towards a smaller frequency as the temperature 
drops below T c , instead of staying intact, as at a nodal 




FIG. 3: Color online. The location of the NSL, defined as 
locus of points for which fermion does not change its energy 
between the normal and the superconducting states. Thick 
(blue) lines are NSL which are detectable by ARPES (the 
peaks above and below T c are located at the same negative 
energy), thin (blue) lines are NSL for which the peak above 
T c is at a positive energy. The hybridization parameter A 
increases from (a) to (d). Dots (red) are nodal points at 
which quasiparticle energy is zero in a superconductor (but 
not in the normal state). The initial evolution of NSL is quite 
involved. We consider it in Appendix B. 



point of a "conventional" superconductor with gap nodes. 
At the same time, one can easily make sure that near the 
larger FS ARPES peak shifts to a larger negative fre- 
quency below T c . As a result, there exist lines in k-space 
along which ARPES intensity does not shift upon cooling 
through T c (in the introduction we termed these lines as 
"no-shift" lines, or NSL). 

In the absence of hybridization, the NSL are radial 
beams along y = a cos 29k = ±1 for which superconduct- 
ing dispersion ui(k) = — (£ 2 + A 2 (l ± y) 2 ) 1 ^ 2 coincides 
with the normal state dispersion — £. Once A becomes 
non-zero, u>(k) splits into two branches, 

4 2 = (A 2 +i/ 2 A 2 +( 2 + A 2 )t2VS, (12) 

where S = £ 2 A 2 + y 2 A 2 (A 2 + A 2 ). We analyzed Q2) 
and found that the NSL rapidly undergo a series of bi- 
furcations (see Appendix B) and for A < A c evolve into 
banana-like loops located in between the two hybridized 
FSs and surrounding the actual nodal points (see Fig. [3]) . 
The NSL persist even above A c , when the actual nodes 
disappear, and should be easily detectable by ARPES (at 
least, in theory). 



B. Non-circular Fermi pockets 

For non-circular electron pockets the quadratic Hamil- 
tonian (J5J) is diagonalized by the angle-dependent trans- 
formation fi = ua + vb, /a = — va + ub, with (ti, v) = 
(cos 0,- sin 0) a nd cos 20 = fo- - 6) 2 + 4A 2 , 
sin2</> = 2X/yJ (£i - £ 2 ) 2 + 4A 2 . The diagonalization 
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FIG. 4: Color online. Hybridization of 2D elliptical FSs. 
Panel (a)- the FSs and the location of nodes (red dots). Solid 
and dashed black lines are the actual, hybridized FSs, and 
the original FSs, respectively. Similarly to the case of circular 
FSs, the pairs of nodal points come closer to each other as A 
increases, and merge and disappear at a critical A c . Panels (b) 
and (c) - the evolution of the NSL with increasing A. Solid 
black lines show the hybridized FSs. Thick and thin solid 
blue lines are NSL detectable and undetectable by ARPES, 
respectively. The 2D analysis is applicable to 1111 materials 
with simple tetragonal lattice structure. For 122 materials 
with body-cetered lattice structure, the situation is somewhat 
different, see Fig. [5] 



yields 



with 



£ 

k 



£,bb\b k 



1/2 



£a,6 = ^(a + &)± [A* + (&-&) 74 
The Nambu matrix M becomes, instead of Eq. 



(13) 



(14) 



£a A(l + a 2 ) Aai ' 

A(l + a 2 ) -f„ Aai 

Aai 6 A(l-a 2 ) 

. Aai A(l-a 2 ) . 



(15) 



where a% = y sin 20, a 2 = y cos 20. The location of the 
nodes is again specified by det M = 0. We have 



dot M = A 2 [(f +&)+aa(6-fa)] 
+ [ Ub + A 2 (y 2 -l)] 2 . 



(16) 



Both terms in (1161) are non-negative and must vanish si- 
multaneously at the locations of the nodes. This sets two 
conditions, one on £ tt) j, and one on the angle k . Solving 
the coupled set we find the same behavior as for circular 
pockets - the eight nodes are located in between the hy- 
bridized electron FSs, and the pairs of neighboring nodes 
move towards each other as A increases, merge along sym- 
metry directions at a critical A c and disappear at larger 
A (Pig. H]). The critical A c is given by 



Xr 



(a 2 -!) 1 / 2 [((£"i-6)/2) 2 + A 



21 V2 



where fi and f 2 are dispersions fi^, &,k along one of 
four symmetric directions 9 k = 0, tt/2, 7r, or 37r/2, and the 
value of k is fixed by the condition fi = £ 2 (a — l)/(a + 1) . 
For elliptical pockets, with fi )2 = fc 2 /2mi j2 + fc 2 /2m 2i i — 

l«, (fi — $2)/2 = aW (1 — a l e l): where e = (mi-m 2 )/(mi + 
m 2 ) is the eccentricity. If the two FSs were nearly circular 
but of different radii before the hybridization (as in body- 
centered tetragonal 122 systems), i.e., £i ^ = fc 2 /(2mi) — 
A*; £,2,k = k 2 / (2m 2 ) — [i, and mi = mil + e), m 2 = m(l — 
e), then fx — £2 — —kpe/(m 2 ), where fc^ ~ \J2m\i. In 
both cases, A c obviously increases when electron pockets 
becomes non-identical. 

We computed the dispersion by solving det(cj(fc) — 
M) = and again found that NSL form banana-shape 
closed loops around the true nodes. We show the results 
in Fig. |4] In similarity to the case of circular pockets, 
NSL survive even when the actual nodes disappear. 



(a) 




(b) ^r-# <^-% <§^©> 



FIG. 5: Color online. Panel (a) - the 3D FSs of 122 systems, 
like BaFe2(Asi- a ,P :!; )2, and the location of vertical loop nodes 
shown by white lines. Two warped cylinders are shown sep- 
arately, but the smaller one is actually inside the larger one, 
as the arrow indicates. Panel (b) - cross section of the actual 
FSs at various k z . Red dots mark the location of the nodes. 



(17) 



IV. THE EFFECT OF HYBRIDIZATION ON 
THE NODES, 3D BODY-CENTERED 
TETRAGONAL LATTICE 

In 122 systems with body-center tetragonal lattice, two 
electron pockets coupled by A are separated by (n, n, 7r), 
where the third component is along k z direction [26| . The 
pairs of pockets at (0,7T, 0) and (w, 0,7r), and at (ir, 0, 0) 
and (0, 7r, 7r) are co-axial ellipses, which are well separated 
even without hybridization [l?} ■ The difference fi — <£ 2 is 
then large, and for k z near and near tt the critical A c 
as given by Eq. (|17l) . is large and generally comparable 
to Fermi energy. On the other hand, at k z = tt/2, the 
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two electron pockets coupled by A are identical, and can 
be well approximated by circles. For these pockets, the 
critical A c is given by (JTTJ) and is quite small, of order A. 
For X c (tt/2) < A < A c (0,7r), the the nodes near k z = tt/2 
are eliminated, and nodal lines form vertical loops, which 
are centered at k z = and k z = tt and close before 
reaching k z — ir/2. We show this gap structure in Fig. [5] 



V. APPLICATION TO P-CONTAINING 
PNICTIDES. 

The three Fe-pnictide materials with hole and electron 
pockets, for which experimental data strongly suggest 
the presence of gap nodes, are LaFeOP, with T c < 5K 
(Ref. 24]), the family BaFe 2 (Asi_ a; P ;c )2 with the high- 
est Tc around 30K (Ref. (H), and LiFeP with T c < 5K 
(Ref. 25]). All three materials contain phosphors. k z - 
integrated probes like penetration depth, thermal con- 
ductivity, specific heat, and NMR [28| all show the be- 
havior consistent with line nodes. In particular, thermal 
conductivity k scales linearly with T at low T and dis- 
plays V~H behavior in a magnetic field, and A(T) — A(0) 
is also linear in T down to very low T. 

The results of ARPES study of t he g ap structure are 
controversial. Laser ARPES study [29| probed the gap 
near the three hole pockets and found all of them al- 
most angle independent, at least for k z probed by laser 
ARPES. These authors argued that the nodes must be 
on electron pockets. Synchrotron ARPES data were, on 
the other hand, interpreted as evidence for a horizontal 
line node at some k z on one of hole FSs. We argue in Ap- 
pendixOthat another fitting procedure of the data from 
Ref. [34 ] . which, we believe, is more appropriate, is con- 
sistent with some k z dispersion of the gap but no nodes 
on hole FSs. And the most recent synchrotron ARPES 
data show no nodes on hole pockets and strong gap 
variation on electron pockets. It appears that more work 
is needed to resolve the structure of the gaps on hole 
and electron FSs in ARPES studies. From theory view- 
point, a horizontal node on a hole FS is possible [26| . 
but less justified than nodes on electron pockets, as the 
latter appear quite naturally due to competition between 
inter-pocket repulsion between hole and electron pockets, 
which favors s" 1 superconductivity, and intra-pocket re- 
pulsion, which is against any superconductivity. The gap 
on electron pockets acquires cos 29 variations to reduce 
the effect of intra-pocket repulsion and allow supercon- 
ductivity to develop 0,13 ■ This reasoning is consistent 
with the argument [6| that a replacement of As by P 
changes the hight of a pnictide with respect to Fe plane, 
which effectively reduces inter-pocket electron-hole inter- 
action, forcing the gap to develop nodes on electron pock- 
ets to reduce the effect of intra-pocket repulsion. We 
therefore believe that nodes more likely reside on elec- 
tron pockets, as we suggest in our analysis. 

The structure of the nodes on electron pockets has 
been discussed in the context of the analysis of ther- 



mal conductivity data in BaFe2(Asi_ x P x )2. Measure- 
ments of the oscillations of thermal conductivity as a 
function of a direction of a magnetic field have been re- 
ported recently [23J, and cos AO component of these os- 
cillations has been interpreted using the same form of 
the gap on electron pockets as in our study: A(fc 2 ) = 
A(l ± a(k z ) cos 28). The best fit to the data yields 
ot{k z ) > 1 for k z near and tt and a(k z ) < I for k z 
near 7r/2. This form of a(k z ) implies that nodes form 
vertical loops centered at k z = and k z = tt. This is 
precisely what we found in our calculations. We there- 
fore argue that our calculation provides microscopic ex- 
planation of the appearance of vertical loop nodes in 
BaFe 2 (Asi_ a; P a: )2. 



VI. CONCLUSION 

To conclude, in this paper we considered how the orig- 
inally nodal s H gap changes if we include into the effec- 
tive tight-binding model for Fe atoms an additional term 
with momentum transfer (tt, tt, 0) in 1111 systems and 
(tt, 7T, tt) in 122 systems. We show that such a term is gen- 
erally present because the hopping between Fe orbitals 
is primarily an indirect hopping via pnictogen (chalco- 
gen) orbitals, and pnictogen (chalcogen) atoms are lo- 
cated above and below the Fe plane in a checkerboard 
order. In 122 systems this order flips between neigh- 
boring Fe planes along z-axis. This additional hopping 
term hybridizes the two electron pockets and affects the 
gap structure. We found that the pairs of neighboring 
nodes (points were the quasiparticle energy is zero below 
T c ) approach each other as hybridization increases and 
disappear once the hybridization parameter A exceeds a 
certain threshold A c . We argued that in 122 systems, like 
BaFe2(Asi_ 2 ,P :E )2, the threshold value depends on k z and 
is much smaller near k z = tt/2 than near k z = and tt. 
In this situation, at intermediate A, the gap nodes form 
vertical loops which are centered at k z = and tt and 
close up before reaching k z — tt/2. 

We also found that in k z cross-sections where the nodes 
are present, they are located away from the hybridized 
FSs. As a result, at a nodal point, the peak in the ARPES 
energy distribution curve shifts, upon cooling through T c , 
from a negative frequency to a smaller frequency. We 
showed that the nodes are surrounded by the "no-shift" 
lines - the subset of k— points at which ARPES peak 
does not shift between T > T c and T < T c . These lines 
initially form beams along the directions where gap 
vanishes, but they rapidly evolve as A increases, and for 
A < A c form banana-shape loops in (k x , k y ) plane around 
the nodes. We propose to search for these "no-shift" lines 
in ARPES measurements. 

We acknowledge useful conversations with D. Basov, L. 
Bascones, S. Borisenko, A. Coldea, D. L. Feng, R. Fer- 
nandes, P. Hirschfeld, I. Eremin, A. Kordyuk, S. Maiti, 
Y. Matsuda, I. Mazin, T. Shibauchi, R. Thomale, M. 
Vavilov, I. Vekhter, A. Vorontsov, and H.H. Wen. We 
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gular dependence of the hybridization amplitude. This 
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(A.C.) 



Appendix A: Microscopic mechanism of 
hybridization 

In this Appendix we present microscopic derivation of 
Eq. and obtain explicit expressions for A(fe) in terms 
of microscopic parameters. We assume, following earlier 
works, that the hopping between Fe atoms occurs via 
pnictogen/chalcogen. Such a process gives rise to two 
types of Fe-Fe hopping terms: the ones with zero mo- 
mentum transfer and the ones with momentum transfer 
Q. The terms with zero momentum transfer (and the 
ones with momentum transfer Q between fermion states 
near the center and the corners of the lFeBZ [ll[) give 
rise to the electronic structure with hole and electron 
pockets. We assume that these terms are already incor- 
porated into the tight-binding Hamiltonian of Eq. ([TJ, 
and focus only on the terms with momentum transfer 
Q, which involve fermions with momenta near (0,7r) or 
(-7T, 0) and give rise to Eq. ([3]). 

Equation ([3]) implies that 



Afc = (l/jl,k\H\lp2,k+Q) 



(Al) 



where H is the full hopping Hamiltonian. The band 
operators ipi,k an d V^.fc+Q are linear combinations of 
atomic orbital operators. We consider d xz , d yz and d xy 
orbitals, and neglect d z 2 and d x i_ y i orbitals which do 
not contribute to the states near the Fermi level Q. Let 
\fs,n) denote orbital operators localized at site n with 
s = 1,2,3 standing for d xz , d yz and d xy respectively. 
The corresponding wave functions are shown in Fig. [5] 
The band operators are 



3 

\^2,k+Q) = ^2l2^k+Q)\fs,k+Q) 



where 



(A2) 



(A3) 



The fc-dependent coefficients 7f(#fe) (7f($/c)) specify or- 
bital contents of the pockets centered at (0,7r) or (ir, 0), 
respectively, and are input parameters for our consider- 
ation. 

Consider momentarily a single layer. We have n = 
(n, m) and Q = (it, it). Using 



J(k+Q)n 





t Px\ 



y 



X 



FIG. 6: Schematic representation of the three Fe orbitals d x 
d yz , d xy and p x and p y orbitals of As/Se. 



we can separate the sums over n into contributions from 
even and odd sublattices, i.e., split J2 n into + 
where the hrst (second) sum is limited to even (odd) 
values of n + m. Equation (|A2[) then takes the form 



\^ k+Q ) = 



E+E 

n n 

I n 

E-E 

n n 



e ifen E^)l/s.n>> 

s=l 

3 

e ikn J2^(0k)\fs, n ). 



(A5) 



Substituting (IA5|) into (|A1|) we obtain 



\k=-Y 



E-E 



n n 



ih{n—n ) 



52[ll(0k)]*Y2(8k)(fs',n>\H\fs 



(A6) 



The sites within each sublattice are identical, hence 
(fs',n'\H\f Sin ) depends only on n' — n. Introducing 

A l s , a = (fs', n +l\H\f s>n ) n + m = 2p 

A'j, s = (f s >,n+i\H\f s , n ) n + m=2p+l, (A7) 

where p is an integer, we re- write (|A6j) as 



A, =~ E 



Al 



7l(0k)}7l(e k )- (A8) 



The summation in Eq. (|A8[) formally runs over all I but 
in reality does not extend beyond second neighbors. The 
same expression is obtained in the 3D case, but then the 
sum over I extends to neighbors in XY plane and along 
z axis. 



(-1) 



(A4) 



Equation (|A8I) expresses in terms of the band struc- 
ture parameters and hopping amplitudes. We see that 
the hybridization parameter is non-zero only if the hop- 
ping amplitudes are different for even and odd sublat- 
tices. 

Further analysis requires the evaluation of the orbital 
hopping amplitudes and the knowledge of the orbital con- 
tent of band operators. As we said, we neglect direct 



hopping between iron atoms and focus on a second or- 
der tunneling processes via pnictogene/chalcogene. The 
amplitudes (|A7j) are then expressed in terms of hopping 
integrals from a d— orbital on a Fe site to one of three 
p— orbitals on a neighboring As/Se site. The p z orbital 
has much smaller overlap integral compared to the two 
in-plane orbitals, and we neglect it. The in-plane or- 
bitals are presented in Fig. [B] Their wave functions are 
maximized in the direction along x[ = (x + y)/V2 and 
x' 2 = (—x + y)/^/2, and we denote then by p x ' ijX ' 2 and 
Px^x'z for pnictogen/chlacogen above and below the iron 
plane, respectively. 

Within the tight-binding approximation the hopping 
parameters between Fe and As/Se sites can be specified 
by a set of overlap integrals 8V a and SV^, as shown in 
Figs. [7] and [5] Here S is the the deviation of As/Se atom 
an Fe plane. The wave functions of d xz and d yz orbitals 
change sign under 8 —5, hence the overlap integrals 
with p Xi and p Xi have opposite signs. 

In what follows we consider first the model with only 
d xz and d yz orbitals. We show that the hybridization pa- 
rameter vanishes in a 2D case and in 3D systems for 1111 
lattice structure, but is generally non-zero for 3D systems 
with 122 lattice structure, except special directions in k 
space where Afc vanishes. We then include into consider- 
ation d xy orbital (i.e., consider three-orbital model) and 
show that in this case the hybridization parameter is non- 
zero already in 2D and in 1111 systems. For 122 systems, 
Afc in the three-orbital model is less anisotropic than in 
the two-orbital model and non-zero along all directions 
offe. 



1. Two-orbital model 

We follow earlier work [l8[ and assume that the pocket 
at (0, 7r) is predominantly made out of d xz orbital, and 
the one at (7r,0) is made out of d yz orbital, i.e., in our 
notations "f{{9k) — S St i and 7f (#fe) = 5 S ,2- 

a. Single Fe layer 
For a single iron layer, Eq. (|A8I) simplifies to 

Afc=5£(4i-4i) e ~ <W - (A9) 
i 

where, we remind, A l 2 x and A 2 i are hopping amplitudes 
from d xz orbital at site n to d yz orbital at cite n + I 
staring from n at even or odd sublattice, respectively. 
Each amplitude describes a two-stage process: the elec- 
tron first hops from d xz orbital to one of the two p or- 
bitals of pnictogen/chalcogen, and then hops from this p 
orbital to d yz orbital on either first or second neighbor, 
see Fig. [§Ia) . Because both d xz and d yz orbitals are odd 
in z, elementary amplitudes of hopping from d to p have 
different sign for hopping originating from even and odd 
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FIG. 7: The overlap integrals between iron orbital d xz , and 
in-plane As/Se p orbitals are denned by two parameters SV n 
and SVcr. We show Fe atom from even sublattice, for which 
As/Se atoms are displaced off the Fe plane by +8 along the 
main diagonal (x — y), and by — <5 along x = —y. The p- 
orbitals at +8 are labeled as p and the ones at —8 as p. We 
assume that 8 <C 1, in which case the overlap integrals are 
linear in 8. (a) overlap between d xz and p x i orbitals, (b) 
overlap between d xz and p x i orbitals. 



sites. However, taken to second order, the amplitudes 
A 2 i and A 21 turn out to be completely equivalent. Wc 
illustrate the equivalence in Fig. [S] for I = (1, 1). As a 
consequence, A\ 1 = A 21 , i.e., vanishes. 

For completeness, we also computed the hybridization 
amplitude between d xz and d yz orbitals for Q = 0: 

Hifl - ]T h h2 (k)fl(k)f 2 (k) . (AlO) 

k 

This term has been computed before @ and is part of the 
bare Hamiltonian which gives rise to hole and electron 
pockets. Using our method, we reproduce earlier result 
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FIG. 8: Same as in Fig. [7] but for d yz orbital of Fe. (a) 
overlap between d yz and p x ^ orbitals, (b) overlap between 
d yz and p x > orbitals. 



for /ii,2, as we now demonstrate. 

Performing the same computations that for Q = (7r, tt) 
lead to Eq. IA91 we obtain for Q = 

^4E(4i+4i)^'- ( A11 ) 
i 

Because A l 21 = A'\ 1 , Eq. (jAllj) becomes 

hi,2 = A,i^ lkl ■ (A12) 
i 

In the second order perturbation theory we have for, e.g., 
1 = (±1,±1) (see Fig. 

,(i,i) _ _ _ ,4(1.-1) - -A ( ~ hl) 

^2,1 — /1 2,1 — ^2,1 _ ^2,1 

= |V CT 2 -^ 2 ), (A13) 




(b) 




FIG. 9: The second order hopping processes within a single 
iron layer shown as a top view (a) and side view (b). The 
equality of the two amplitudes, A™' 1 ' = A'* 1 ' 1 ' follows from 
the observation that both orbitals, d xz and d yz axe odd under 
the reflection z — > —z. 



where the energy denominator E g is the energy separa- 
tion between 3d to 4p orbitals. Substituting Eq. (|A13|) 
into Eq. (jXT2]) we obtain 

h h2 = -^{Va - V*) sinfc, sink y + . . . , (A14) 

where dots stand for the contributions from other I. This 
result is in ful l ag reement with earlier calculations, see 
e.g., Refs. HUM 



b. 3D crystals, 1111 materials 

For 3D systems with 1111 structure, Fe layers at dif- 
ferent z are al equivalent, and 3D folding vector Q = 
{■n, 7r, 0). The arguments displayed in previous section ap- 
ply to this case as well, i.e., for two orbital model = 0, 
even if we include into consideration inter-layer tunnel- 
ing. 
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c. 3D crystals, 122 materials. 

In 122 materials with body-centred tetragonal crys- 
tal structure the situation is qualitatively different. The 
first observation is that the hybridization vector is Q — 
(tt, tt, tt) because even and odd sublattices are formed by 
Fe atoms located at n = (n, m,p) with n + m + p even 
or odd, respectively. This can be also be understood by 
noticing that the Fe-only lattice is simple cubic, but, be- 
cause of As/Se, Fe lattice has an foe structure with the 
basis, or alternatively a rock salt structure. The folding 
vector Q — (tt, tt, tt) appears as additional Bragg peak 
due to the transition from a simple cubic to an fee lat- 
tice. Without interlayer tunneling, A& is still zero, but 
inter-layer tunneling makes it finite as we show below. 

Inter-layer tunneling in real systems is a complex pro- 
cess which at least partly involves Ba atoms (Ref. (171]). 
We will avoid this complication and consider a toy model 
in which there is a direct tunneling between pnicto- 
gen/chalcogen atoms located at the same (x, y) in differ- 
ent layers. Because the position of pnictogen/chalcogen 
atoms relative to Fe plane oscillates along z direction, 
there are two different tunneling amplitudes for such 
processes: the one between pnictogen/chalcogen located 
above n-th plane and below (n + l)-th plane, and the 
other between pnictogen/chalcogen located below n-th 
plane and above (n + l)-th plane (see Fig. [T0)l . 

Another peculiarity of 122 systems is that the eccen- 
tricity of electron FSs in 122 systems changes sign be- 
tween k z = and k z = tt. i.e., the FS at (0,7r0) is 
elongated along the same direction as the FS at (tt, 0, tt). 
In the full model, this is due to the change of the rela- 
tive weight of d xy orbital between k z = and k z = tt 
(Ref. [17J). In two-orbital approximation, we model the 
change of sign of eccentricity by requiring that d xz and 
d yz orbitals interchange between k z = and k z = tt [i.e., 
we require that at k z = 0, the pocket at (0, tt) is made out 
of d xz orbital and is elongated along y-axis, and the one 
at (tt, 0) is made out of d yz orbital and is elongated along 
x— axis, while at k z = tt, the pocket at (0, tt) is made out 
of d yz orbital and is elongated along x, and the one at 
(tt, 0) is made out of d xz orbital and is elongated along 
y\. The hybridization term with Q = (tt, tt, tt) then con- 
nects (0, tt) pocket at k z = and (tt, 0) pocket at k z = tt, 
which have the same orbital character, i.e., s and s' in 
the Eq. (|A8I) for A& are the same. 

The process which gives rise to a non-zero A& is a three- 
stage process in which an electron from d xz or d yz orbital 
on an even sublattice in layer n hops to one of two p- 
orbitals on a pnictogen/chalcogen, then hops vertically 
to a pnictogen/chalcogen located near next Fe layer, and 
then hops to the same d— orbital on that layer. One can 
easily verify that the terms in A& in which do not can- 
cel out between even and odd sublattices are the ones 
with I = (±1,±1,±1). There are two different hopping 
amplitudes for these momenta (see Fig. [T0|) : for even 

sublattice we have Aji' 1 ' 1 -' = Ais' -1 ' ^ = Als = 
_ , ,(-1,1,1) _ ,(1,-1,1) _ ,(1,1,-1) _ 

si-ss — ^M? Jiss — ^iss — s±ss — 




FIG. 10: Interlayer tunneling processes leading to a finite 
hybridization, (|A15[l in 122 material. 



Ai s l ' 1 ' = A-2, where s — 1, 2, and for odd sublattice 
,'(l,i,i) _ ,'(l-i,-l) _ .'(-i,-i,i) _ '(-i.i-i) _ 

Slss — ^*ss — -tt-ss — slss — si2, 

/(-M,i) _ ,'(1,-1,1) _ ,'(1,1,-1) _ ,'(-1,-1,-1) _ , 

-fiss — Ass — Siss — jiss — A-l • 

Substituting these amplitudes into Eq. IA81 we see that 
in the real part of A^ the contributions from even and 
odd sublattices cancel out, but in the imaginary part of 
Afe they add up, such that 

Afc = — 4i(A% — A%) sink x smk y smk z . (A15) 

The expression (IA15[) is odd in all three momenta. This 
result is expected because the three reflection symmetries 
have been broken by hybridization. The overall factor 
of i reflects the fact that time reversal symmetry is not 
broken (i.e., A fe = Al fe ). 

For completeness, we also computed the hybridization 
term between d xz and d yz orbitals and found it is also 
non-zero, but this time A& is real and is an even function 
of fe. 



2. Three orbital model 

We next consider how A& changes if we consider more 
realistic situation when d xy orbital also contributes to 
the states near the FS. 



a. 2D case 

We use as an input the results of previous calcula- 
tions [l|, which found that the FS at (0, tt) is constructed 
out of d xz and d xy orbitals, and the one at (tt, 0) is 
constructed out of d yz and d xy orbitals. To a reason- 
able approximation, 7x(#fe) = sin#, 7i(0fe) = cos#, and 
7| (6) = 7{(7r/2 — 8), where 9 is counted from y-axis (as 
written, the formulas are valid in the first quadrant, for 
< 6 < tt/2). Accordingly 

ipi,h = fi.k sin 6 + fz.k cos 6 
$2,k+Q = h.k+Q cos 6 + / 3 ,fc+Q sin 6 (A16) 

The hybridization amplitude between /i.^ and /2,/e+Q 
vanishes, as we found before, but we show in this sec- 
tion that the amplitudes between f± t k an d fs,k+Q and 
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between /^fc and f2,k+Q are non-zero. In our notations, 
we then have 



Afc = Ay + A 



3,1 , x 2,3 
k ' 



where 



Af=sin 2 ^(4, 1 -4i) 



-ikl 



and 



Af = cos 2 0^(4,3-^ 



(A17) 



(A18) 



(A19) 



Because the wave function for d xy orbital is even under 
the reflection z — > —z an the ones for d xz and d yz orbitals 
are odd, A 31 — —A l 31 and A' 2 l 3 = —A 23 . Equation 
(|A19|) then reduces to 



A fc = \(sin 2 9 A l 3 l + cos 2 9A l 23 )e 



i 



(A20) 



We now need to prove that this expression is non-zero. 
As before, we consider hopping between Fe sites as a 
two-stage process via neighboring p-orbitals. This pro- 
cess gives rise to hopping to nearest and next nearest 
neighbors on the iron lattice, i.e., we need to consider 
I = (Ixi ly) with l x ,l y — 0,±1. The overlap integrals 
with the p orbitals on As/Se have been defined in Figs. [7] 
and [5] In explicit form 

(Px' 1 ,n+(l/2)(x+y)\H\fl. n ) = —6V a , 
(Px' 2 ,n+(l/2)(x+y)\H\fl,n) = SV„, 
(px' 1 ,n+(l/2)(x+y)\H\fl,n) = SV a , 
(Px' 2 ,n+(l/2)(x+y)\H\fl,n) = SV K , (A21) 

where, we remind, 5 is the deviation of As/Se from an Fe 
plane. 

We define overlap integrals between d xy and p orbitals 
in Fig. QT] as 

(f3,n+x\H\Px 2 ,n+(l/2)(3s+§)) = G, 
(f3,n+x\H\Px' lt n+(l/2)(3s-§)) = G, 
{h,n+x\H\p x ' 2 ,n+(l/2)( X +y)) = G, 
\f3,n+si\H\px' 1 ,n+{l/2)(x-y)) = G (A22) 

Other overlap integrals between d xy and in-plane 
p— orbitals vanish, 

Using these notations, we obtain after a simple algebra, 



^3\i 0) =(fl,n\H\Px' 3 ,n+(l 



(l/2)(x+y)> 



X (Px' 2 ,n+(l/2)(x+y)\H\f 3 ,n+x) 



-{fl,n\ H \Px' 1 ,n+(l/2)(x-y)) 
x -^-(Px[,n+(l/2)(x-y)\H\f3, n +x) ■ (A23) 






FIG. 11: Overlaps of iron d xy orbital with As/Se orbitals p x ^ 

(a) and p x i^ (b). In the case (a) the tunneling is possible only 
along the main diagonal direction, x — y, while in the case 

(b) it occurs only along the perpendicular direction x — —y. 



A *3,l^ - (h,n\H\p x ' ltn +(i/2)( 



x+y)l 

X -^r(Px' 1 ,n+(l/2)(x+y)\H\.f3,n+x+y) ■ (A24) 

and similar results for other A 3 1' ly ^ and A 23 ■ Sub- 
stituting the expressions for the overlap integrals, Eqs. 
(|MT|) . CMU), we obtain 



.4 



(±i,o) _ 2SVttG „(o,±i 



3,1 



E„ 



) ^3,1 — u > 



(A25) 



i(±i,±i) _ ^SV a G 4 (±i,=f1) _ ^SV a G 



A?~'=T^, A^' +1 '=T^, (A26) 



and similarly 

,(0,: 
1 2,3 



,(o,±i) _ 28V V G (±i i0 ) _ 

^i' i iTi i ^2,3 — 



0. 



(A27) 
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,(±i,±i) _ | 8V a G 



A 



(±l,Tl) 
2,3 



T§5. (A28) 



Subs titut ing Eqs. ([A25]) . (|A26)) . (|A27j) and (|A28j) in 
Eq. (|A20p and summing over four nearest and four next- 
nearest neighbors, we obtain 

Afc = ^— [V^r (sin 2 6* sin fc^ + cos 2 sin k y ) 



E, 



+14- (sin 9 sin k x cos fc y — cos 9 sin fc y cos k x )] . (A29) 



Introducing small k = k - 
a small-size pocket cos 9 



(0,7r), using the fact that for 
= k y /\k\, sin9 = k x /\k\, and 



sin &x j sin kf. 



-k y , and extending the analysis to 



other quadrants (i,.e., to negative k x and k y ), we finally 
obtain 



At- = 



AiSG 



(V, - V„) (\k~\l - \k\l) . (A30) 



We see that hybridization parameter A& is generally non- 
zero, except for diagonal directions k x = ±k v . This 
agrees with the result of numerical calculations [19] . The 
authors of Ref . Qj^ found that in the presence of spin-orbit 
interaction |A| remains non-zero even along the diagonal 
directions. 



b. 3D case 

In 1111 systems the hybridization vector is Q = 
(tt, 7r, 0), and for weak inter-layer tunneling Afc weakly 
depends on k z and is nearly the same as in Eq. IA29I 

In 122 systems the hybridization vector is Q = 
(tt, tt, tt), and the full expression for A& in 3D case is the 
sum of d xz — d yz contribution, Eq. (|A15[) . and the contri- 
bution from the processes involving d xy , Eq. (IA29[) . The 
two-orbital contribution is strongly fc z -dependent, while 
the three-orbital contribution is fcz-independent if we ne- 
glect inter-layer hopping, and weakly depends on k z if 
we include it. Because the two contributions vanish along 
different symmetry directions, the total hybridization pa- 
rameter is non- vanishing for all k, though it has minima 
at k z — 0, tt, and along k x — ±k y and k x — 0, k y = 
in a given cross-section of k z (k x , k y are measured with 
respect to a center of a pocket). The anisotropy of A& is 
further reduced if we add spin-orbit interaction (l7l[T9j . 



Appendix B: Evolution of the nodal lines at small 
hybridization 

In this Appendix we discuss the initial evolution of the 
NSL from radial beams directed transverse to the FSs, to 
tangential lines, directed along the FSs. We present the 
results for circular pockets. The evolution of the NSL for 
elliptical pockets is quite similar. 

In the normal state, the dispersions for hybridized cir- 
cular pockets are = ± A, where = (k 2 — 
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FIG. 12: Color online. The initial evolution of the NSL de- 
fined as locus of points for which fermion does not change its 
energy between the normal and the superconducting states. 
Dashed (black) lines show the FS before the hybridization. 
Thick (blue) lines are NSL which are detectable by ARPES 
(the peaks above and below T c are located at the same neg- 
ative energy), thin (blue) lines are "nodal lines" for which 
the peak above T c is at a positive energy. The hybridization 
parameter A increases from (a) to (f). Dots (red) are nodal 
points at which quasiparticle energy is zero in a superconduc- 
tor (but not in the normal state). Panel (a) A < A c ,i, panel 
(b) A = A c ,i, panels (c,d) A c ,i < A < A c ,2, panel (e) A = A c ,2, 
panel (f) A > \ c ,2- Both A c ,i and \ c ,2 scale as A 2 /fi and are 
much smaller than A c ~ A, at which the nodes disappear. 
The evolution at A < A c is shown in Fig. [31 



kp)/(2m) is the dispersion in the absence of hybridiza- 
tion. In the superconducting state the dispersion is given 
by Eq. (TIT 



s (Jfe) = (A 2 + y 2 A 2 + £ 2 + A 2 ) T 2v^ , 



(Bl) 



where S = £ 2 A 2 + y 2 A 2 (A 2 + A 2 ) and y = acos29 k . At 
A = 0, Eq. (|B1[) reduces to a conventional expression 
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uj\ 2 (k) = £ 2 + A 2 (l ± y) 2 . The dispersions in the normal 
and superconducting state coincide along 4 radial lines 
specified by y = ±1. Along these directions, s"* 1 gap has 
accidental nodes on one or the other electron FS. 

Once A becomes non-zero, the lines transform into 
banana- type loops still elongated transverse to the FS, 
and "domes" at k > kp (Fig. IT2"a ) The loops close up 
at k = 0, and at k = kp and y « ±1. The NSL then 
undergo several topological changes at A ~ A 2 //i, which 
is much smaller than critical A c ~ A. At the first critical 
A = A c i the bananas touch pairwise along the directions 
9 k = 0, ±7r/2 and tt (Fig. At larger A, eight ba- 

nanas transform into four configurations, which resemble 
rabbit ears (Figs. [T2lc-d). The value of A c ,i is 



A c ,i 

At the next critical 



1 



A^jl a 2 + 1 



A 2 
4 M 



(B2) 



(B3) 



the "rabbit ears" touch each other along another set of 
symmetry directions, 9k = ±7r/4,±37r/4 (Fig. [T2"e). and 
at A > A Cj2 detach from k = ((Fig.[12f). The evolution 
of the NSL at larger A ~ A is discussed in the main text 
and is shown in Fig. |4] 



Appendix C: Fitting of the ARPES data 

In this Appendix we discuss the fitting procedure 
used in Ref. [34[ to extract the gap structure from 
ARPES data. The authors of this work symmetrized 
the measured photoemission intensity to get rid of Fermi 
function and extract the spectral function A(k, uj) = 
(—1/tt)G (fci?,w), and used the convention that the gap 
is finite when A(k,uj) at k — kp has two resolved max- 
ima at finite frequencies, and zero when the maximum 
in the spectral function is at zero frequency. The same 
procedure was used earlier to identify Fermi "arcs" in 
the cuprates (36[. This fitting procedure requires care 
because the maximum in A(kp,uj) can be at zero fre- 
quency even when the gap is non-zero. This happens 
when fermionic damping is finite (as it always is at a 
finite T, even in a conventional s-wave superconductor) 
and is larger than the gap. 

To extract the gap from the data, the authors of [34[ 
related the spectral function to fermionic self-energy in a 
standard way, as 



A(k F ,oj) 



1 



E (k F ,Lj) 



7T ((J - Z'(kp,Uj) 2 + (£"(fcf,w)) 2 

and modeled the self-energy at each k z as 

A 2 



T,(k F , w) — —ij + 



UJ 



(CI) 



(C2) 




FIG. 13: Color online. The spectral function A{kF,uj) (a.u.), 
Eq. (|C1|) . with self-energy given by (a) Eq. (|C2|) and (b) 
Eq. (EU). The solid black line is for 7 = 0.5 A. The dashed 
blue line is for 7 = 3. OA. In the case (a), the spectral func- 
tion has a dip at oj — and a peak at a finite frequency for 
any non-zero A. In contrast, in the case (b), the peak in 
the spectral function can be at uj — despite that the gap is 
non-zero. This happens when the broadening is large enough, 
7 > \/3A. Notice that in this situation the peak at uj — is 
very broad. 



If one uses this fit, one finds that the maximum in the 
spectral function is at uj = only when the superconduct- 
ing gap vanishes (see Fig [T3h). Indeed, the self-energy 
given by (|C2|) diverges at uj — and therefore the spectral 
density has a dip at uj = for any nonzero A. Because 
the measured spectral function a hole pocket is peaked 
at uj = near k z = tt, the authors of [34| concluded that 
the gap on the a pocket must vanish for this k z . 

We argue that the fermionic self-energy of a dirty BCS 
superconductor contains the damping term 27 not only 
as the stand-alone constant but also in the denominator 
of the A 2 term, along with uj, i.e., the true self-energy is 



T, t (k F ,uj) = -i'y + 



uj + 47 



(C3) 



This form of the self-energy has been extensively dis- 
cussed in the context of the physics of Fermi arcs in the 
cuprates [37l - [39| . Using this self-energy, one obtains that 
the spectral function has a maximum at uj = even when 
the gap is finite, provided that 7 > \/3A. We show in 
Fig. IT3"b the spectral function obtained using the self- 
energy from Eq. (|C3I) with a non-zero A. We see that 
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the maximum in A{kF,uS) is either at a finite frequency 
or at u) — 0, depending on the interplay between A and 7. 
Observe that, when the maximum is at to — 0, it is quite 
broad, much like in the experimental data in [34| near 



k z = 7r. We argue therefore that the data of the ARPES 
study in [34j are in fact consistent with anisotropic but 
still no-nodal gap along the a pocket. 
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